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Abstract -The generation of entanglement between two oscillators that interact via a common 
reservoir is theoretically studied. The reservoir is modeled by a one-dimensional harmonic crystal 
initially in thermal equilibrium. Starting from a separable state, the oscillators can become entan- 
gled after a transient time, that is of the order of the thermalization time scale. This behaviour is 
observed at finite temperature even when the oscillators are at a distance significantly larger than 
the crystal's interparticle spacing. The underlying physical mechanisms can be explained by the 
dynamical properties of the collective variables of the two oscillators which may decouple from or 
be squeezed by the reservoir. Our predictions can be tested with an ion chain in a linear Paul 
trap. 
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' Introduction. — Dissipation and decoherence seri- 
ously limit the possibility to observe quantum effects in 
the macroscopic world. Understanding their origin and 
dynamics at the microscopic level is a fundamental issue 
fbr quantum technologies, which aim at achieving control 
of the dynamics of scalable quantum systems. Microscopic 
models based on Hamiltonian dynamics have been stud- 
ipd over the last century. A textbook example is Brownian 
motion [Hll]: damping and thermalization of a particle's 
ip-Otion is here found when the particle is embedded in a 
crystalline structure, which acts as a thermal bath under 
specific conditions [5H1]. The physical picture obtained 
from these models supports the understanding that loss 
of coherence emerges from the creation of quantum corre- 
lations between system and reservoir [SHJ- 

Recent works pointed out that the interaction with a 
common external reservoir may endorse the creation of 
entanglement between two physical systems, for instance 
two spins [HHH] or two oscillators P^HT5] which arc not 
directly coupled to each other. These studies show that 
the bath can exert an active role in establishing entangle- 
ment. In most cases, this observation can be explained 
using symmetry reasons: a collective variable of the sys- 
tems is decoupled from the reservoir such that the corre- 



sponding eigenstates form a Decoherence-Free Subspace 
(DPS) [10]. Provided that the systems initial states and 
the reservoir temperature fulfil certain conditions, entan- 
glement can then be found at time scales at which a sin- 
gle system would have otherwise thermalized. A different 
mechanism for entanglement generation is based on the 
direct coupling of a collective variable with the bath [19j 
and results from the squeezing of its variance at sufficiently 
low temperatures of the reservoir • 

The dependence of entanglement when the systems are 
coupled to distant positions of the reservoir has recently 
been object of several studies [TUMTSlfrTlfTB] . In the ther- 
modynamic limit, entanglement is found between a pair 
of oscillators of a one-dimensional chain only when these 
are nearest-neighbours [131 HZ]. In ref. [TH], it was ar- 
gued that entanglement between two oscillators which cou- 
ple to a bath described by a generalised Caldeira-Leggett 
model [H] disappears at distances larger than the wave- 
length associated with the cutoff frequency of an Ohmic 
spectral density [5] , which corresponds to the interparticle 
spacing when the reservoir is a crystalline structure. Prom 
these works one may expect that this length constitutes a 
limit, above which entanglement cannot be established by 
a reservoir of arbitrarily large size. 
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In this Letter we show that these expectations are to 
a large extent unfounded. We analyse the creation of en- 
tanglement between two harmonic oscillators that interact 
via a common reservoir — a chain of harmonic oscillators 
with nearest-neighbour coupling. Wc first examine the 
conditions for which the chain acts as an Ohmic bath, in 
the setup of fig. [T^) which is analogous to Rubin's model 
when a single oscillator couples to the chain 013]. In this 
particular case, our results are in qualitative agreement 
with ref. [19| . where the bath was described in terms of 
an Ohmic spectral density. We then demonstrate that 
the same reservoir can also support the creation of en- 
tanglement at long times between oscillators which couple 
to different, distant particles of the chain, as depicted in 

fig. Eh). 

Hamiltonian dynamics and initial states. Our 

model consists of two oscillators with mass frequency 
ri, position and momentum (a = 1,2) that in- 
teract with a particle of a linear chain. The linear chain 
is composed by 2N particles of mass m, position Xj and 
momentum pj (j = ±1, . . . , ±N) that couple via nearest- 
neighbour interaction. The edge particles are pinned by 
a harmonic trap with frequency ujb- The total Hamilto- 
nian H = Hs + Hb + Hj comprises the Hamiltonians of 
the two oscillators, the chain, and their mutual interac- 
tion: 
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Here, n and 7 describe the coupling strengths and 
OJi = (jJb {^i,-N + <^i,Af) sets appropriate boundary condi- 
tions on the chain. The Heaviside function Q{t) indicates 
that Hi vanishes before the time t = 0. For later conve- 
nience, we introduce the cutoff frequency Wcut = a/ 4k/to 
of the chain's normal modes and denote by a the interpar- 
ticle distance, with 2a being the wavelength of the mode 
at frequency Wcut |3]. Moreover, we define a reference fre- 
quency VLq such that = (1 -I- e) JIq with a dimensionless 
parameter e that will be specified later. 

The dynamics of the composite system is evaluated nu- 
merically, assuming that each oscillator is initially pre- 
pared in a squeezed vacuum state of Hamiltonian ([T]) . We 
denote by r the squeezing parameter which is taken to 
be real such that the variances are AX^(O) — c^'^'^a^ /2 
and AP2(0) = e2''ft7(2a2) with a = ^h/{Ivmo). The 
chain is prepared in a thermal state at temperature T, 
and its density matrix reads ps(0) = jZ with 

Z = Tr{e-^^«} and ^ = (kBT)-^. By integrat- 
ing the equations of motion of the composite sys- 
tem, we determine the time-dependent covariance ma- 

trix ^;. = i(ce+ec)-(C)(C) (a*,;^^ {1,2, 3, 4}) 
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Fig. 1: Microscopic models a) and b) 
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in terms of the dimensionless variables of the os- 
cillators {X[,P{,X^,P^) with X'^ = X„/a and 
= Per a/h. The entanglement between the oscilla- 
tors is evaluated by means of the logarithmic negativity 
En = max{0, — ln(2i'_)}, where £>- is the smallest sym- 
plectic eigenvalue of the partially transposed covariance 
matrix V = AV^'A that results from a mirror refiection 
A = diag(l, 1, 1, -1) with respect to P^ j^ . 

Microscopic model of an Ohmic bath. Wc first 
focus on the model depicted in fig. [1^) where both oscilla- 
tors couple to one edge of the chain, such that j = k = N. 
A simple transformation shows that only the centre-of- 
mass (COM) motion of the two oscillators with coordi- 
nate X+ = (Ai +X2)/\/2 couples with particle N, while 
the relative motion with coordinate A_ = (Ai — X2)/V2 
is decoupled. The coupling of the COM to the chain is 
analogous to Rubin's model of a single defect in a crys- 
tal jl| . Thermalization of the defect is expected when the 
crystal acts as an Ohmic bath. In order to recover this be- 
haviour, we study the spectral density of the chain [HHI] 



N _2 



2m ^ — ' (jjj 



(3) 



where 7^ is the strength of the coupling between the 
COM oscillator and the chain's normal mode at frequency 
uji, whereby the normal modes diagonalize the quadratic 
Hamiltonian Hb + jx'j^ ■ 

By choosing the masses and the parameters k, 7 and 
u!b according to the prescriptions discussed in [J] , we find 
J+{i-o) ~ Lu for io < n, which corresponds to the spec- 
tral density of an Ohmic bath and is displayed in fig. [5^). 
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The thcrmalization process is numerically investigated by 
analysing the behaviour of the covariance matrix for times 
smaller than the revival time ircv ~ L/cs, where L — 2Na 
and Cs ~ Wcuta/2 are the length and the sound velocity 
of the chain, respectively. Note that trcv grows linearly 
with N. For 2N = 2500 and 7 = ClMf^g, we observe 
that all correlations between COM and relative coordi- 
nates vanish at sufficiently long times t (but t < irov), 
and the only non-vanishing components of the covariance 
matrix are the constant widths AX^, AP^, and the el- 
ements AXl, API and i (X_P_ -|-P_X_). The latter 
oscillate at the frequency 2rt.y where fi^ = 77 is the fre- 
quency of the oscillator for the relative motion which en- 
tails the frequency shift due to the interaction Hamilto- 
nian Hj through the parameter 77 = y^l+^y/jAKW) . In 
accordance, the quantity fAr(i) — — ln(2i>_(t)) oscillates 
with 20.^ between the values f]^"'^ < S'^''''. 

Depending on the initial squeezing parameter r and the 
temperature T, we encounter three different situations for 
the entanglement of the two oscillators when the COM mo- 
tion has reached its "stationary" state. Our findings qual- 
itatively agree with the results of ref. [TH], where this be- 
haviour was predicted by means of a non-Markovian mas- 
ter equation with an Ohmic (sub-Ohmic, super-Ohmic) 
spectral density [23]. The different entanglement be- 
haviours are summarised in fig. [2}d) where the logarith- 
mic negativity is displayed as a function of r and T. 
They are indicated using the nomenclature introduced 
in [19) : (i) entanglement sudden-death (SD) [24] when any 
transient entanglement disappears, such that 8^^^ < 0, 
(ii) entanglement sudden death and revival (SDR) when 
< < S^'^^, and (iii) entanglement no sudden 
death (NSD) when 5™'" > 0. 
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Fig. 2: a) Spectral density J+{u}) for the model in fig. [T^), 
and b) contour plot of the logarithmic negativity EN{r,T) at 
long times (but smaller than the revival time) as a function 
of the squeezing parameter r and temperature T (in units of 
/if2o/fcs). The black region (NSD) indicates the parameter 
regime in which the oscillators are entangled (the SDR phase 
is not visible in this graph, but located between SD and NSD). 
The parameters are 2iV = 2500, e = 0, m = M/2, lob = ^2^0, 
K = Mnl, 7 = O.lMf^o and ojcut = 2V^Qo. 

Entanglement between the two oscillators is found when 
the quadratures of the COM and relative motion exhibit 



two-mode squeezing [5S]. An analysis of the logarithmic 
negativity at long times shows that the oscillators are en- 
tangled when either inequality (I) AXf (T) < ^e^l''"''^! 
or (II) APf (T) < §e-2|'-'-sl is satisfied, providing a rela- 
tion between T and r. Here we also introduced the squeez- 
ing parameter^ = hirj. Two-mode squeezing has two 
different origins. The first one is connected to the DPS for 
the relative motion which leads to entanglement for large 
\r\ according to inequality (I). At sufficiently low temper- 
atures and for r « rg, however, is squeezed due to 
the coupling with the bath [2]|7] and entanglement is found 
when (II) is fulfilled. The corresponding region is located 
where the two NSD-regimcs meet in fig. ^jp) (not visible 
for the chosen scale). 

Long-distance entanglement. We now turn to 
the generation of entanglement when the oscillators are 
coupled to different particles of the reservoir, as illus- 
trated in fig. [T)d). Let the two particles of the chain be 
separated by a distance d that is larger than the inter- 
particle separation a. We set j = —n and fc = ?i in 
eq. such that d = {2n + l)a. After performing the 
transformation xf = {xi ± X-i)/ V2 (i = I, . . . , N) for the 
bath particles, the whole dynamics reduces to the time 
evolution of two separate harmonic oscillators with coor- 
dinates X± that interact with two independent reservoirs. 
Figure [3^) displays the spectral density of the reservoir 
coupling with the centre-of-mass {J+(u;)) and the relative 
coordinate (J_(w)). One observes that both spectral den- 
sities monotonically increase for frequencies to sufficiently 
close to zero. For oscillator frequencies il that lie in this 
interval, the reservoir acts as an Ohmic environment for 
each collective variable. In this case, the oscillators reach 
a stationary state which exhibits no entanglement. This 
behaviour is found at any distance d > a for sufficiently 
small values of 51 and 7, and is in agreement with the re- 
sult of ref. [18] , where the reservoir was described in terms 
of an Ohmic spectral density. 

A significantly different behaviour is obtained when the 
(shifted) oscillator frequency Vl-y takes on values, for which 
the spectral density oscillates. In particular, when we 
choose e such that coincides with a frequency value 
Lo ~ uj^ at which one of the spectral densities vanished, 
J±i^E) =0, the corresponding variable X+ or X^ de- 
couples from the bath. In this case, entanglement at long 
times exhibits similar dynamics as for the model in fig.[T^), 
but with an important difference: the decoupled mode in- 
volves both the variable X± as well as the variables xf of 
the particles interposed between the oscillators. 

Figure [SJa) displays the logarithmic negativity at long 
times for d = 9a, showing that long-distance entanglement 
occurs for sufficiently low temperatures T. Remarkably, 

^An initially squeezed state with squeezing parameter rg corre- 
sponds to the ground state of an oscillator with trap frequency Qj. 
^The zeros of J±{i^) read Ldjj, i_ = (jJcut sint/)^ with 

</>± = If (2fc± - i T I), fc± e {1, . . . , ± ij} and [x] as floor 
function. 
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Fig. 3: a) Spectral densities J+{uj) (dashed line) and J-(w) 
(solid line) corresponding to the centre-of-mass and relative 
motion of the oscillators, respectively, for the model of fig. [TJj) 
with d — 9a and b) contour plot of the logarithmic negativity 
at long times as a function of r and T. Here, e = —0.086, 
2N = 1500; all other parameters are as in fig. [21 



we observe that the two oscillators cannot be entangled 
above a certain temperature T, independent of the ini- 
tial squeezing \r\. This is in sharp contrast to the action 
of the Ohmic environment (see fig. Hb)): where entangle- 
ment can be generated at any temperature, provided |r| 
is sufficiently large. 

The reason for this difference can be traced back to the 
type of DFS which corresponds to each model. In order to 
clarify it, let us recall the physical picture which explains 
entanglement in the model of fig. la). In this case, only 
the COM coordinate X+ of the two oscillators interacts 
with the reservoir, whereas the relative coordinate X- is 
decoupled. FigureH^) provides a sketch of the correspond- 
ing model in terms of collective variables. In this situation 
the squeezing of the initial state of the relative motion is 
preserved (even though the squeezed quadrature rotates 
at the oscillator frequency), while the correlations with 
the COM vanish. The relative motion will hence exhibit 
squeezing at long times, provided that the two oscillators 
are initially prepared in two independent squeezed states. 
This will lead to entanglement (two-mode squeezing) at 
any temperature as long as the initial squeezing parameter 
r is sufficiently large. In the model of fig. lb), in contrast, 
both the COM and relative coordinates of the two oscil- 
lators couple to a reservoir defined by the coordinates xf 
and x~ , respectively, with xf = (xt ± X-i)/V^- However, 
in this case the two reservoirs are independent one from 
the other. This situation is illustrated in fig. HJd). Here, 
one can identify a mode, which involves (X_ ) and the 
corresponding bath particles with coordinates Xi,... x^ 
{x^ , . . .x~), and that decouples from the rest of the chain 
provided that the spectral density J+{u!) {J-{uj)) vanishes 
at the shifted oscillator frequency fl^. The initial state of 
this mode results from the initial squeezing of the corre- 
sponding collective variable of the two oscillators as well 
as from the thermal state of the involved chain oscillators. 
As a result, the width of the Gaussian state associated 



with the oscillator X+ (AT-) of the decoupled mode de- 
pends also on the initial temperature of the chain. This 
property leads to the disapperance of long-distance entan- 
glement above a certain temperature, irrespectively of the 
initial squeezing. 
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Fig. 4: a) Microscopic model of fig.[T^) and b) of fig.HJj) here re- 
ported in COM and relative coordinates X± = [Xi ± X-2)l\f2 
of the two oscillators. In a), whereas the COM motion couples 
to the harmonic chain, the relative coordinate decouples from 
the rest of the dynamics. In b) COM and relative motion of 
the oscillators couple to two different reservoirs defined by the 
coordinates xf = (xt ± X-i)/^/2. 

We further remark that our numerical simulations show 
that entanglement decreases linearly as the distance d is 
increased. At fixed distances d, however, the logarithmic 
negativity remains constant when the number N of par- 
ticles in the chain is increased. Hence, entanglement at 
large distances can b found at arbitrarily long times in 
the thermodynamic limit. 

Long-distance entanglement crucially relies on the non- 
Markovian nature of the bath and it appears when the 
trap frequency of the system particles is tuned to values 
at which the spectral density vanishes. Nonetheless, the 
behaviour of an Ohmic bath can be re-established by re- 
moving one of the system particles in fig. [TJd), in which 
case the remaining oscillator thermalizes with the har- 
monic chain. This feature already highlights an impor- 
tant difference between our microscopic model and sys- 
tems that are commonly employed for long-distance cor- 
relations such as dipoles that couple via a photonic-band 
gap material [lTJ[26]. 

We note that the choice of our model is based on Ru- 
bin's harmonic crystal [3] , which serves as a paradigm for 
the analysis of quantum Brownian motion [5]. Unlike the 
systems studied in refs. [13l[T7] and [18], where no dis- 
tant entanglement is found, our model (i) possesses no 
discrete translational invariance and (ii) is characterised 
by a spectral density which is Ohmic only for a small fre- 
quency interval, whose width decreases as the oscillator's 
distance is increased. This shows in particular that ap- 
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proximating the environment with an Ohmic (sub-Ohmic, 
super-Ohmic) spectral density does not catch the rele- 
vant symmetries of the Hamiltonian which can support 
a DFS and thus entanglement creation. Deviations from 
this ideal behaviour, such as spatial inhomogencity, are 
therefore expected to be detrimental as they destroy the 
symmetries of the coupling. 

Our prediction could be tested with an ion string in a 
linear Paul trap, building on the proposals in [57] ■ The 
oscillators could be identified with the radial vibrations 
of two heavier ions embedded in the chain, provided that 
they are sufficiently distant such that their mutual repul- 
sion can be neglected. The axial and radial degrees of 
freedom of a defect correspond in our model to one oscil- 
lator and one particle of the chain, respectively. They can 
be harmonically coupled by means of dipolar potentials. 
With adequate radial trap frequencies and mass ratios be- 
tween the defects and the other ions, the radial oscillators 
can be decoupled from the radial excitations of the rest of 
the chain. Moreover, their frequency could be tuned with 
respect to the axial spectrum by choosing appropriately 
the trap aspect ratio [28] . 

Summary. — We have studied a paradigmatic model 
based on Rubin's harmonic crystal which supports entan- 
glement generation between two harmonic oscillators that 
indirectly couple via a one-dimensional harmonic chain. 
The creation of long-distance entanglement presented in 
this letter is neither a result of a controlled dynamical 
evolution, nor is it due to the finite size of the reservoir as 
in ref. fl5| . Entanglement is found in the thermodynamic 
limit — > oo and at finite temperatures (i) for distances 
between the oscillators that are significantly larger 
than the interparticlc spacing of the chain and (ii) for 
times at which a single oscillator would have otherwise 
thermalized. These results may find useful applications 
in quantum technologies, for instance within the context 
of hybrid quantum networks based on quantum systems 
at the interface to a solid state environment. 
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